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Abstract 

ABSTRACT: We investigate the mechanism of spin gap formation in a 
two-dimensional model relevant to Mott insulators such as CaV^g. From 
the perturbation expansion and quantum Monte Carlo calculations, the origin 
of the spin gap is ascribed to the four-site plaquette singlet in contrast to the 
dimer gap established in the generalized dimerized Heisenberg model. 
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In previous papers, [0-0] we have investigated the mechanism of spin gap formation 
in a class of two-dimensional generalized antiferromagnetic Heisenberg (AFH) models with 
dimerization. From this analysis, the origin of the spin gap in low-dimensional systems such 
as the spin-Peierls system, the ladder system and the Haldane system has been identified as 
the dimer gap in a unified way. 

Recently, Taniguchi et al. found that the temperature dependence of the uniform mag- 
netic susceptibility and the nuclear magnetic relaxation time T\ of CaV/tOg indicate the 
existence of the spin gap@ In this paper, we investigate the possible mechanism of the spin 
gap in CaV/tOg. 

The crystal structure of VO5 pyramid layer in CaX^Og is shown in Fig. |l](a). Since the 
valence of V atom is 4+, V atom can be treated as the localized spin with S = 1/2. In the 
case of copper oxide systems such as high-T c materials and Sr n _ 1 Cu„ + i0 2n , the electrons 
are localized in the c? 7 orbitals on Cu 2+ sites. However in this system, the Fermi level lies 
in the d e orbital and hence the localized spins should mainly come from electrons in the 
d e orbitals on V 4+ sites. The degree of the orbital degeneracy with three-fold degenerate 
d e orbitals has not been clarified yet, because the anisotropic crystal field due to possible 
Jahn- Teller distortion has not been estimated precisely. In any case, the antiferromagnetic 
exchange coupling is expected between the adjacent d e orbitals on the nearest-neighbor V 
atoms. We also expect that the exchange coupling between the nearest neighbors is dominant 
due to the d e character. Based on these facts, the S = 1/2 AFH model with the network 
shown in Fig. |](b) is a good starting point to discuss the Mott insulating phase. The 
Hamiltonian is written as 

H = JY,Si-Sj, (1) 

where denotes the nearest-neighbor bonds. 

We first calculate the spin gap by quantum Monte Carlo (QMC) simulation to check 
whether our model is relevant or not. Here the spin gap A s is defined as the energy difference 
between the lowest state with S* ot!il = J2i Si = 1 and the singlet ground state. The algorithm 
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of the QMC simulation is based on the world-line method combined with the Suzuki- Trotter 
method || and the checker-board decomposition. The temperature T scaled by J is 
taken as 0.1, which is sufficiently low to observe the ground state properties for the present 
purpose. The Trotter numbers N T are 40, 60, 80 and 100. The number of QMC steps is 
(1 ~ 5) x 10 7 for each Trotter number. The extrapolation formula used to obtain A s from 
the finite N T is given by A S (L) = A s (oo) + A/N 2 . The lattices we have calculated have 
sizes lxl, 2x2, 3x3 and 4x4 cells under the periodic boundary condition, where a 
cell includes 20 sites. The value of the spin gap extrapolated to the thermodynamic limit is 
estimated to be (0.11 ± 0.03) J with the form of the fitting function A s + A/L 2 , as is shown 
in Fig. |[ 

We also calculate the temperature dependence of the uniform susceptibility for 2 x 2 cells, 
which is shown in Fig. |[ The susceptibility in the thermodynamic limit should be close to 
this result, since the the spin gap for 2x2 cells is close to the value in the thermodynamic 
limit. As a reference, the susceptibility of the conventional square-lattice AFH model on 
8x8 lattice is also shown in Fig. |3], which may describe the thermodynamic properties 
within statistical errors. [[J This also shows the existence of the spin gap in our model. 

Next we investigate the origin of the spin gap. The four-site plaquette is treated as 
the unperturbed Hamiltonian while the other bonds are taken as the perturbation. As is 
shown in Fig. |], the strength of the spin exchange couplings in the four-site plaquette is 
taken as J, while the strength of the others is J'. The ground state of the unperturbed 
Hamiltonian is the product state of singlets on the four-site plaquettes, in other words, a 
resonating valence bond state in each plaquette. The first excited states consist of the triplet 
state on one of the plaquettes and the singlet on the others. The degeneracy is lifted by 
the first-order perturbation expansion (PE) through the hopping of the triplet due to the 
translational symmetry. Then the energy of the triplet state has a wave-number dependence. 
We calculate the energy of the ground state and triplet states using the second-order PE. 

(2) 

Then the ground state energy per site e G and the energy difference between the ground 

(2) 

state and the triplet state A s (fc) are obtained as 
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1 47 J 72 

A s 2) (fc) = J + - J' (cos k x a + cos k y a) - — — 
6 504 J 

1 J' 2 1 J' 2 

+— — (cos 2 k x a + cos 2 k y a) — -— cos k x a cos k y a, (3) 
54 J a J 

where a is the lattice constant between plaquettes, as is shown in Fig. The dispersion of 
the triplet excitation A^\k) is shown in Fig. [|. The lowest triplet excitation is located at 
k = (it, 7r). An interesting point is that the spin gap survives (0.205 J) even if the value of J' 
is equal to J. On the other hand, if we take the unperturbed Hamiltonian as J' bonds with 
J I J' taken as the perturbation, the ground state energy per site given by the second-order 
PE is obtained as 

4 2) = -^'d + ~). (4) 

('2') f(*2,} 

When J' is equal to J, e G = — 0.537J obtained from eq. (g) is lower than e c; = —0.469 J 
in eq. The dispersion of the triplet state given by the first-order PE in terms of J/ J' is 
described as 

Ag 1 -* (k) — J' — - J(cos k x a + cos k y a) . (5) 

The perturbation breaks down and the spin gap closes at J/J' = l even in the first-order 
PE. 

The above results imply that the origin of the spin gap is basically the four-site plaquette 
singlet state rather than the dimer singlet. The spin gap in CaX^Og has been estimated 
from the susceptibility as A s ~ 100K, whereas the gap amplitude obtained from the 
QMC is A s /J ~ 0.11. Although at present we do not have available data for J in CaA^Og, 
one may argue that the calculated value of the spin gap appears to be smaller than the 
observed value because the exchange interaction in vanadium oxides is expected to be at 
most several hundred K. One possible origin of the gap enhancement is the frustration effect 
arising from the next-nearest-neighbor exchange coupling. It is also noted that the effect of 



the orbital degeneracy and the orbital correlation effect could be important for a quantitative 
estimation of the gap. 

In order to investigate the mechanism of the spin gap in more detail from the viewpoint 
of the four-site plaquette singlet formation, we investigate the one-dimensional analog of 
this model. The lattice structure is shown in Fig. |6](a). We calculate the spin gap by 
the exact diagonalization (ED) method. The size dependence of the spin gap is shown in 
Fig. |7|. After the extrapolation to the thermodynamic limit, we estimate the spin gap as 
A s ~ 0.60J, which is larger than the value of the spin gap in the ladder model. |lO|Jl~T|| The 



fitting function is the same as the one in the previous analysis. In the second-order PE, the 
triplet state energy is obtained as 

1 31 J' 2 1 J' 2 

A s (k) = J+ -J'cosk x a - + \Qg~ cos2kxa - ^ 

When J' is equal to J, the value of the spin gap is A s = 0.658J at k — 7r, which is close 
to the estimation from the ED. These results in the one-dimensional model show that the 
four-site plaquette singlet is a good starting point for discussing the ground state in a class 
of lattices constructed from the four-site plaquettes connected by a small number of bonds. 
This one-dimensional model may be not a toy model but a relevant model in some transition 
metal oxide compounds if a lattice structure such as that in Fig. |6](b) is realized. 

In summary, we have studied the mechanism of the spin gap in CaA^Og. From the QMC 
as well as the PE, the origin of the spin gap formation is found to be the four-site plaquette 
singlets. 
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FIGURES 



FIG. 1. (a)The lattice structure of VO5 pyramid layer in CaV^g. Open circles represent V 
sites, while solid circles show oxygen sites. (b)The network of the S = 1/2 AFH model. In the 
QMC calculation, we take 20 sites as a cell as shown by the wavy-line square. 

FIG. 2. Size dependence of the spin gap in our model. 

FIG. 3. Temperature dependence of the uniform magnetic susceptibility for 2 x 2 cells. Solid 
circles represent the susceptibility of our model for 2 x 2 cells, while open circles show the suscep- 
tibility for the conventional square lattice AFH model for 8 x 8 lattices. 

FIG. 4. The lattice structure in the perturbation calculation. Bold lines represent the bonds 
with J and thin lines those with J' . The distance between the nearest-neighbor four-site plaquettes 
is a. 

FIG. 5. The wave-number dependence of the triplet excitation energy in the second-order 
perturbation calculation. The inset shows the route of the wave vector taken for the abscissa. 

FIG. 6. (a)The one-dimensional four-site plaquette model. The definitions of the sym- 
bols are the same as in Fig. 1(b). (b)The lattice structure of a possible quasi-one-dimensional 
four-site-plaquette system of transition metal oxides. 

FIG. 7. Size dependence of the spin gap in the one-dimensional model. The sizes we have 
calculated are 8,16 and 24 sites. 
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